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The stochastic (^*-theory in d— dimensions dynamically develops domain wall structures within 
which the order parameter is not continuous. We develop a statistical theory for the y*-theory 
driven with a random forcing which is white in time and Gaussian-correlated in space. A master 
equation is derived for the probability density function (PDF) of the order parameter, when the 
forcing correlation length is much smaller than the system size, but much larger than the typical 
width of the domain walls. Moreover, exact expressions for the one-point PDF and all the moments 
{ifi") are given. We then investigate the intermittency issue in the strong coupling limit, and derive 
the tail of the PDF of the increments tp[x2) — tp{xi). The scaling laws for the structure functions of 
the increments are obtained through numerical simulations. It is shown that the moments of field 
increments defined by, Cb = {\'p{x2) — ip{x\)\^), behave as \x\ — X2\^'' , where ^i, = & for 6 < 1, and 
^6 = 1 for 6 > 1. 
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I. INTRODUCTION 

Thirty five years ago Wilson and Fisher [1] emphasized 
the relevance of the (y9^-thcory to understanding the crit- 
ical phenomena. Since then, the theory has become one 
of the most appealing theoretical tools for studying the 
critical phenomena in a wide variety of systems in statis- 
tical physics. In the strong coupling limit, the (/3^-theory 
develops domain walls, a phenomenon which is of great 
interest in the classical and quantum field theories [2-10]. 
The dynamical (^^— theory - what is usually referred to 
as the time-dependent Ginzburg-Landau (GL) theory - 
provides a phenomenological approach to, and plays an 
important role in, understanding dynamical phase transi- 
tions and calculating the associated dynamical exponent 
[11-16]. The time-dependent GL theory for superconduc- 
tors was presented phenomenologically only in 1968 by 
Schmid [17] (and derived from microscopic theory shortly 
thereafter [18]), when the first modulational theory was 
derived in the context of Rayleigh-Benard convection 
[19,20]. Moreover, the GL equation with an additional 
noise term has been studied intensively as a model of 
phase transitions in equilibrium systems; see, for exam- 
pie, [21]. 

In the present paper we consider the stochastic (/s*- 
theory in the strong coupling limit. This limit is singular 
in the sense that, the equation that describes the dynam- 
ics of the system develops singularities. Therefore, start- 
ing with a smooth initial condition, the domain-wall sin- 



gularities are dynamically developed after a finite time. 
At the singular points the field ^p{x, t) is not continu- 
ous. We derive master equations for the joint probability 
density functions (PDF) of ip and its increments in d di- 
mensions. It is shown that in the stationary state, where 
the singularities are fully developed, the relaxation term 
in the strong coupling limit leads to an unclosed term in 
the PDF equations. 

Using the boundary layer method [22,23], wc show that 
the unclosed term makes no finite contribution (anomaly) 
in the strong coupling limit, and derive the PDF of and 
its moments, ((/s"), in the same limit. Wc also investigate 
the scaling behavior of the moments of the field's incre- 
ments defined by, 5lp = ip{x2) — ip{xi), and show that 
when \x2 — xi\ is small, fluctuations of the ip field have 
a bi-fractal structure and arc intermittent. The inter- 
mittency implies that the structure function defined by, 
Cb = {\(p{x2) — (p{xi)\'') , scales as la;2 — a; i ]'''', where Ch is a 
nonlinear function of b. It is also shown numerically that 
the moments of the field's increments, {\ip{x2) ~ (p{xi)\^) , 
behave as, \x2 ~ a^il^*", where Xh = b for b < 1, and ^6 = 1 
for 6 > 1. 

The rest of this paper is organized as follows. In the 
next section we present the model that we wish to study, 
and describe some of its properties by solving it numer- 
ically. In III and IV we derive master equations for the 
order parameter of the model, and for the field's incre- 
ments and its PDF tail. The numerical simulations for 
extracting the scaling exponents are described in V. The 



paper is summarized in VI, while the Appendices provide 
some technical details of the work that we present in the 
main part of the paper. 



II. THE MODEL AND THE COUPLING 
CONSTANT 

The standard GL (p^-theory describes a second-order 
phase transition in any system with a one-component 
order parameter Lp{x) and the ip ^ ~ip symmetry in 
a zero external field. The theory is described by the 
following action, 
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(1) 



where r = T — T^, with being the critical temperature, 
and k is the diffusion coefficient. We consider the case 
in which r > 0. For d > 2, the critical temperature 
Tc is finite, while in one dimension (ID), Tc = 0. The 
parameter g characterizes the strength of the fluctuation 
interaction, or the coupling constant. The equation of 
motion is given by. 



fcV Lp + Ttf - 



6^ 



0. 



(2) 



The critical dynamics of the system is described by a 
stochastic equation of a particular form - the Langevin 
equation (for a comprehensive review on Langevin type 
equations see [24] )- given by 
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ry(x,t) 



(3) 



where 77(x, t) is a Gaussian-distributed noise with zero 
average and the correlation function. 



(7?(x, i)7y(x', t')) = DqD{^ - ^)6{t - t') , 



(4) 



with £'(x— x') being an arbitrary smooth function. Typi- 
cally, the spatial correlation of the forcing term is consid- 
ered to be a delta function in order to mimic short-range 
correlations. Here, though, the spatial correlation is de- 
fined by 



D(x 



TTcr- 



■ exp 



(x-x') 



/^2 



(5) 



where a <^ L endows a short-range character to the ran- 
dom forcing. It is useful to rescale Eq. (3) by writ- 
ing, ip' ~ f/f^, x' = x/xo, and t' = t/to. If we let 
to = 1/r, ifo = (6r/5)l/^ and, xq = [idog)/i6T^)]^/'', 
all the parameters are eliminated in Eq. (3) except for, 
k' = [6/{Dog)]'^^'^T'^/''-^k, and one finds that 



where k' is now the effective coupling constant of the the- 
ory. The weak and strong coupling limits of the theory 
are then defined, respectively, by, /c' — > oo and fc' — > 0. 
In the weak coupling limit one can use numerical simula- 
tions and the Feynman diagrams to calculate the critical 
exponents. On the other hand, to solve the problem in 
the strong coupling limit we need other techniques to de- 
rive the stochastic properties of the fluctuation field [16]. 
The nonlinearity of Eq. (6) in the strong coupling limit 
gives rise to the possibility of formation of singularity in 
a finite time. This means that there is a competition 
between the smoothing effect of diffusion (the Laplacian 
term) and the (p^ term. Let us now describe the main 
properties of the GL theory in the limit, k' — > 0. 

i) The unforced GL model [r]{x,t) ~ 0], with given 
initial conditions, develops singularities in any spatial di- 
mension. In one spatial dimension (ID) the singularities 
are developed in a finite time tc as k' —>■ 0. At such singu- 
lar points the field (p, representing an order parameter, is 
not continuous. In 2D the unforced GL model develops 
domain walls, characterized by singular lines with finite 
lengths (that depend on the initial condition). Under 
these conditions, the field (p is discontinuous when cross- 
ing the singular lines. In three and higher dimensions the 
structure of the singularities can be more complex. For 
example, in 3D the singularities are domain walls where 
the field ip is discontinuous. 

In Figs. 1 we show the time evolution of the or- 
der parameter p of the unforced GL model in 2D, in 
the hmit k' [Eq. (6) with 77 = 0]. We have 
used the finite-element method to numerically solve the 
Langevin equation with k' —^ and the initial condition, 
(p{x,y,0) = sin X sin?/. Such initial conditions are typi- 
cal, and were used only for simplicity. The time scale for 
reaching the singularity is of the order of fc'~^/^. As Figs, 
la and lb indicate, it is evident that at times t < tc (in 
the limit, k' 0) the (p field is continuous. At t = t,. the 
(p field becomes singular; see Fig. Ic. 

(ii) Similarly, for a forcing term which is white noise 
in time and smooth in space, singularities are developed 
in any spatial dimension in the strong coupling limit and 
in a finite time, t,, = tc.r;, as k' —^ 0. For example, in 2D 
the boundaries of the domain walls are smooth curves. 
In Fig. 2 we demonstrate the time evolution of the order 
parameter ip of the forced GL model in 2D in the limit 
k' —!■ 0. Starting from a smooth initial condition, as 
shown in Figs. 2a and 2b, it is evident that for times 
t < tc^n the p field is continuous. At t = tc,i7 the field 
becomes singular; see Fig. 2c. 



III. MASTER EQUATION OF THE ORDER 
PARAMETER 



dip 
'dt 



k'\/'^p + p-ip^ + ij{x,t) , 



(6) 



In this section we derive a master equation to describe 
the time evolution of the PDF P{p, t) of the order pa- 





FIG. 1. Snapshots of the order parameter at various 
times in the 2D unforced (/3*-theory in the strong coupling 
limit, k' — » 0. Top two figures show the ip field before the 
singularity develops. In the bottom figure, which is for time 
scales greater than the time at which singularity develops, 
the ip{x, y) field is not continuous. The initial condition is, 
fix, y, 0) = sin x sin y. 

rameter (p. Defining a one-point generating function by, 
Z(X) = (0), where 6 is defined by, Q = cxp[—iXip{x,t)]. 
Using Eq. (6), the time evolution of Z is governed by 

Zt = -i\k'(V'^ifexp[-i\tp{x,t)]) - iX{ipQ) 

+iX{f^Q) - i\{r]<d) - \^k{{))Z , (7) 

where, k{x) = D()D{x.), and we have invoked Novikov's 
theorem (see Appendix I), which is expressed via the re- 
lation, 

{t] cxp[~iXif{x, t)]) = -iXk{0)Z . (8) 

Now, using the identities, —iX{ipcxj>[—iXip{x,t)]) = XZ\, 
and —iX{ip^ cxp[—iX(p{ic, t)]) — XZ\XX, the generating 
function Z satisfies the following unclosed master equa- 
tion 
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FIG. 2. Snapshots of the order parameter </p of a random- 
ly-driven 2D (p^'-theory in the strong coupling limit. In the 
simulations the relation between the forcing length scale a 
and the sample size L is, o" ~ L/3. The forcing strength Do 
is 0.1. 

Zt = -iAfc'(VVexp[-a(y9(x,t)]) + XZx + XZxw 

-X^k{{))Z . (9) 

The -iXk'{V'^ipex-p[-iXip{x,t)]) term of Eq. (9) is the 
only one which is not closed with respect to Z. The 
PDF of order parameter P{lp) is constructed by Fourier 
transforming the generating function Z: 

Pi^,t)^ J ^cxp(*A^)Z(A,0 . (10) 

Thus, 

Pt = -[i^-^^)P]^ + k{0)P^^ 

— Xexp{iXip){V'^ipexp[—iXip{x,t)]) . (11) 
2tt 

It is evident that the governing equation for P{ip, t) is 
also not closed. 



Let us now use the boundary layer technique to prove 
that the unclosed term [the last term of Eq. (11)] makes, 
in the strong coupling limit, no contribution to the gov- 
erning equation for the PDF [22,23]. We consider two 
different time scales in the limit, k' 0. (i) Early stages 
before developing the singularities {t < tc^,^), and (ii) 
in the regime of established stationary state with fully- 
developed sharp singularities (t > tc,ri)- 

In regime (i), ignoring the relaxation term in the gov- 
erning equation for the PDF, one finds, in the limit 
k' 0, the exact equation for the time evolution of the 
PDF for the order parameter (see below for more de- 
tails). In contrast, the limit fc' ^ is singular in regime 
(ii), leading to an unclosed term (the relaxation term) 
in the equation for the PDF. However, we show that the 
unclosed term scales as fc'-'^/^, implying that this term, in 
the strong coupling limit, makes no finite contribution or 
anomaly to the solution of Eq. (11). It is known for such 
time scales (the stationary state) that the (/?-field, which 
satisfies the Langevin equation, gives rise to discontin- 
uous solutions in the limit, k' 0. Consequently, for 
finite a the singular solutions form a set of points where 
the domain walls are located, and are continuously con- 
nected. We should note that k'(pxx, in the limit k' 0, 
is zero at those points at which there is no singularity. 
Therefore, in the limit k' only small intervals around 
the walls contribute to the integral in Eq. (11). Within 
these intervals, a boundary layer analysis may be used 
for obtaining accurate approximation of ip(x,t). 

Generally speaking, the boundary layer analysis deals 
with problems in which the perturbations are operative 
over very narrow regions, across which the dependent 
variables undergo very rapid changes. The narrow re- 
gions, usually referred to as the domain walls, frequently 
adjoin the boundaries of the domain of interest, due 
to the fact that a small parameter (fc' in the present 
problem) multiplies the highest derivative. A powerful 
method for treating the boundary layer problems is the 
method of matched asymptotic expansions. The basic 
idea underlying this method is that, an approximate so- 
lution to a given problem is sought, not as a single ex- 
pansion in terms of a single scale, but as two or more 
separate expansions in terms of two or more scales, each 
of which is valid in some part of the domain. The scales 
are selected such that the expansion as a whole covers 
the entire domain of interest, and the domains of valid- 
ity of neighboring expansions overlap. In order to handle 
the rapid variations in the domain walls' layers, we de- 
fine a suitable magnified or stretched scale and expand 
the functions in terms of it in the domain walls' regions. 

For this purpose, we split tp into a sum of inner solu- 
tion near the domain walls and an outer solution away 
from the singularity lines, and use systematic matched 
asymptotics to construct a uniform approximation of if. 
For the outer solution, we look for an approximation in 
the form of a series in k\ 



y, = <^°"t^^o + fcVi+0(fc'2) 
where ipo satisfies the following equation 



(12) 



(13) 



Indeed, (po satisfies Eq. (6) with k' = 0. Far from 
the singular points or lines, the PDF of (po satisfies the 
Fokker-Planck equation, with the drift and diffusion co- 
efficients being, D^^^ipo,t) ^ tpQ — ipl, and D^^^ = k{0), 
respectively. Reference [16] gives the solution of the time- 
dependent Fokker-Planck equation with such drift and 
diffusion coefficients. At long times and in the area far 
from the singular points or lines, the PDF of (po will have 
two maxima at ±1. This means that we are dealing with 
the smooth areas in Fig. 2c in the stationary state. 

In order to deal with the inner solution around the 
domain walls, we consider the x component normal to 
the domain wall or singularity line, and decompose the 
operator as dxx + ^d-i- the strong coupling limit, 
k' — !■ 0, the term V^_j^(/3 makes no contribution to the 
PDF equation, whereas the term dxx'P' is singular. To 
derive the long-time solution of Eq. (6), we rescale x to 
z = ^^^= and suppose that complete solution of Eq. (6) 
has the form, (p{z, t) = f{z, t) + ta.nh{z). All the effects of 
the initial condition and time-dependence of fix, t) will 
then be contained in f{z,t). We now rewrite Eq. (6) 
with the new variables to obtain 



dtf{z,t) = ld. 



f{z,t) + f{z,t)-f{z,t) 

-3f{z, t) tanh(z) [f{z, t) + tanh(z)] 

+ y^i-l{z,t) . (14) 

The last term of Eq. (14) is zero in the limit k' — > 0. 
Multiplying Eq. (14) by f{z,t) and integrating over z, 
one finds that, 

dt [dzf{z,t) = - f dz{dJ{z,t)f + 2 [dzf{z,t) 



(15) 



where a{z) ~ tanh(z). We show in Fig. (3) the time vari- 
ations of / dzp{z, t) verses t with different types of initial 
conditions. The results show that J dzp{z,t) vanishes 
at long times. Therefore, Lp{z,t) — s- tanh(z), in the limit 
of a stationary state. Let us now compute the contribu- 
tion of the unclosed term in Eq. (11) in the stationary 
state, i.e., 

-k' /"^zAe'^^(vVe^*^^^"'*^) 
27r 



-2 j dzf{z,t)~6 J dzf{z,t)a{z) 
-6 / dza^{z)f{z,tf , 



= -k' 



dX 
2^ 



k'{V'fS[ip-ip{x,t)]), 



(16) 



In the second line of Eq. (16), we have replaced iX with 
differentiation with respect to ip and in the third line the 
integration of A has been carried through. Now. assuming 
ergodicity, the term k' {'V'^(p5((p~ip{x, t))) is converted to, 

= —k' lim — ( dxdvd-i^'^^p5[ip — ip{x,t)\ . (17) 

y— oo V Jy 

In the limit fc' — * 0, only at points where we have singu- 
larity the above term is not zero. Therefore, we approach 
the domain walls' regions as 

= -k' lim ^y^f dxdvd-i{fxx + ^1-iV) 

xJ[(^-(^(x,i)] , (18) 

where fij is the space close to the domain walls. There- 
fore, Eq. (16) is written as 

= -k' lim -^V" / dxdvd-Wxx^\Hi ~ •^{x,t)\ . (19) 

Changing the variables from a; to 2 and integrating 
over dvd-i, we have 

= lim V / edz\(p^J[(p - (p{z,t)] 

= -— lim V I dz(p^J[ip - (p{z, t)] 



where e = {2k'y^^. Assuming statistical homogeneity, 
one finds, 

k' N X Vd 1 r^°^ 
= lim — — !- / dzip^^Slip" Lp{z,t)] (20) 

where N is number of singular lines. The quantity 
NVd-i/V is the density of the singular lines, and fc'/e = 
(fc72)i/2. In the limit, V 00, we denote the density of 
the singularities by p. Therefore, 

= -(^72)1/^/ dz^,,5[p~^{z,t)]. (21) 



Now, by changing the integration variable form z to 
we can determine the integral exactly. We find that, 



-t-00 



+1 



dzip:^z5['f - '^{z,t)]^ / (5[(^ - (^(z,t)] . 

J-l 

(22) 

Using Eq. (6) in the limit, i ^ 00, we determine tpzz 
and (fz. in terms of (f. Multiplying Eq. (6) by (fz and 
integrating over z, we obtain. 



(23) 



where C is an integration constant. In the limit, z — > 
±cx), Lpz = ^ = il- Therefore, C ~ 1/2, and ipzz/'fiz is 
written as, 



w 

1) 



~ ^Fiip X signiip^ — 1) . 
The integral in Eq. (22) is now given by, 

dip'^^5[Lp — (p{z, t)] 
-1 Pz 
+1 

dip{z)\/2ip{z) X sign{ip{z)'^ 

-1 

2ip sign{Lp^ - 1)6'(1 - ip^). 



l)5{^~^{z,t)) 



Now, the unclosed term in Eq. (16) is written as 

-k'{V'^ip5[ip^ip{x,t)])^ = 
= -{k'/2Y'^p{A + 2^/29,2^(^2 _ ^)|^ 

where, A = 2 V26'(l - ip^)e{Lp^ - 1) - V26'(l - cp^). There- 
fore, in the limit k' 0, the master equation takes on 
the following form 

= -[((^ - Lp^)P]^ + k{Q)P^^~ 

{k'l2fl^p\A + 2\/2(^2^(^2 _ ^ (24) 

where we set Pt = in the stationary state. The PDF 
P is continuous at ip = ±1, but its derivative is not. 
By integrating Eq. (24) in the interval [1 — e, 1 -I- e] (or 
[—1 — e, — 1 + e]), one finds 



{k'fl'^p 

m 



(25) 



In the limit, fc' — > the derivative of the PDF will also be 
continuous. Considering the factor (fc')^''^ (24), we 

conclude that, in the strong coupling limit, the unclosed 
term is identically zero. This means that there is no 
anomaly or finite term in the strong coupling limit in the 
master equation for the PDF of the order parameter ip. 
The stationary solution of Eq. (24), in the limit, k' 0, 
takes on the following expression. 



P.i 



N exp 



21 



4fc(0) 



where the normalization constant is given by 



1 

N 



1 



V2 



exp 



1 



8fc(0) 



Ki 



1 



8fc(0) 



(26) 



(27) 



K/r 




10--° I , , „ , 

10^ 10-' 10"- ^ 10' 10" 10' 



FIG. 3. Time variations of J dzf^{z,t) vs. t with dif- 
ferent types of initial conditions. The resuhs show that 
J dzf^{z,t) ^ at long times. 

where the Ka(h) is the modified Bessel functions. To 
derive the moments of ((^") in the stationary state, we 
multiply Eq. (25) by (/s" and integrate the result over <p 
to obtain 

+?^(n-l)A:(0)n((/5"-2^) ^0. (28) 

Equation (28) is a recursive equation for computing all 
the moments in terms of the second-order one. Direct 
calculation then shows that, 

(sfeW) 

and all the odd moments vanish, ((^^'^+^) = 0. Therefore, 
using Eqs. (28) and (29), we are able to derive all the 
moments of the order parameter in the d— dimensional 
Ginzburg - Landau theory in the strong coupling limit. 

IV. MASTER EQUATION FOR THE 
INCREMENTS AND THEIR PDF TAIL 

In this section we derive the PDF and the scaling 
properties of the moments of the increments, P{(p(x2) — 
(p{xi)), for the (/3*-thcory in the strong-coupling limit. 
Defining the two-point generating hmction by, Z(X) ~ 
(0), where Q is defined as 

the time evolution of Z is related to that of (p by 

Zt - -^Al((pt(xl,^)e-'^l^(^l^*)-*^^'^(^^■*)) 
-~iX2{ipt{x2,t)e-'^"^^'='''^-'^''^'-'-'^'*^) . (31) 

Substituting (pt(xi,t) and ipt{x2,t) from Eq. (6), the 
governing equation for the generating function satisfies. 



Zt — AiZai + AiZajAiAi + ^2Z\2 + A2.^A2A2A2 

-{Xl + Xl)k{0)Z -2XiX2k{x)Z , (32) 

Where, x ~ xi ~ X2, and wc have used the fact that in 
the strong-coupling limit, the Laplacian term makes no 
contribution to the PDF equation (see Appendix II for 
more details). Moreover, we have invoked the generalized 
Novikov's theorem for the two-point generating function, 
according to which [23] (see also Appendix I), 

-ai(r;(a;i)e-'^i'^(^i^*)-''^^'^(^^'*') 
-jA2(77(a;2)e-'^i^(^i^*)-'^^'^(^^'*)) 
= -{Xl + Xl)k{0)Z -2XiX2k{x)Z . (33) 

Fourier transforming Eq. (32), the governing equation of 
the joint PDF will be given by 

Pt{ipi,^2) = - V'l)PU 

+2k{x)P^,^^, . (34) 

It is useful to change the variables as, ipi = [w — u)/2, 
and, ip2 = (w + u)/2, and, therefore, d/dipi = d/dw — 
d/du, and, dj dip2 ~ d/dw + d/du. Now, Eq. (34) can be 
written as 

Pt{w,u)^-{wP)^-{uP)u + 

[^{w^ + 3wu^)P]^ + [^{u^ + 3uw^)P]u 

2fc(0) (-P^uu, 4- Puu) + 2fc(a;) (i-'mto Puu ) • (35) 

To derive the governing equation for the PDF of the 
increments, u = Lp2 — ^pi, we integrate over w to find that, 

Pt{u) = -{uP)u + J(u'P)„ 

^{u{w^ I u)P)u + 2[A:(0) - k{x)]Puu , (36) 

where we have used the fact that the joint PDF P{'w, u) 
can be written as, P(w \ u)P{u). It is evident that we 
cannot derive a closed equation for the PDF of m. Indeed, 
to determine P{u) we need to know the conditional av- 
eraging {uP' I u). However, one can derive the tail (both 
the left and right ones) of the PDF in the limit, u — > oo. 
To determine the tail we note that only near the sin- 
gularities, in small separation in space, one finds a large 
difference in the field and, hence, large u. On the other 
hand, near such points or lines, the field w will be very 
small. Therefore, in the limit u — > oo, we can ignore the 
conditional averaging to find that, 

lim (w^ I u) ~ . (37) 

It — *oo 

Therefore, in the limit, u s- oo, wc obtain the following 
behavior for the tails of the P{u) in the stationary state. 



P,,(,,_,,^„,.„p{-J|lM_}. (38) 

To derive the scaling behavior of the moments (w") one 
needs to know the entire range of the behavior of the 
increments' PDF. Here, we are able to only derive the 
equation for the shape of the PDF tails. In the next sec- 
tion, we investigate by numerical simulation the scaling 
behavior of the moments (w") vs. the separation x. 

V. SCALING EXPONENTS OF THE MOMENTS: 
NUMERICAL SIMULATION 

To calculate numerically the scaling behavior of the 
moments with x, when x « 1, we shall use here the 
initial-value problem for the two-dimensional Langcvin 
equation, Eq. (6), in the limit, k' — > 0, when the force is 
concentrated at discrete times [17-20]: 

/(x,y,t) = ^/,(x,2/) (39) 

3 

where both the "impulses" fj{x,y) and the "kicking 
times" tj are prescribed (deterministic or random). The 
kicking times are ordered and form a finite or an infinite 
sequence. The impulses are always taken to be smooth 
and acting only at large scales. The precise meaning that 
we ascribe to the dynamical Langevin equation with such 
a forcing is that, at time tj^ the solution 1^9(2;, y, t) changes 
discontinuously by the amount fj{x,y), 

'Pix,y,tj+) = (p{x,y,tj^) + fj{x,y) , (40) 

whereas between tj^ and the solution evolves ac- 

cording to the unforced ip'^ equation, 

dtip = (p-ip^ + k'V^ip. (41) 

Without loss of generality, we may assume that the 
earliest kicking time is, tj^ = provided that we set, 
fjo = /o: and, ip{x,y,tjg_) = ip{x,y,tQ_) for t < to. 
Therefore, starting from tg, according to Eq. (40) we 
obtain 

'p{x,,y,tQ+) = ip{x,y,ta^) + fo{x,y) , (42) 

and beyond that up to ti^, according to Eq. (41), 

y, h-) = (1 + h)ip{x, y, to+) - (p^{x, y, to+), 
h^h-to (43) 

where, h = ti — Iq. 

It is clear that any force f{x,y,t) which is continu- 
ously acting in time can be approximated in this way by 
selecting the kicking times sufficiently close. Hereafter, 
we shall consider exclusively the case where the kicking is 
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FIG. 4. Exact solution of the PDF of (/j^'-theory and its 
comparison with numerical calculation. 

periodic in both space and time. Specifically, we assume 
that the force in the equation is given by 

+00 

f{x,y,t)=g{x,y) ^ S{t-jT), (44) 

j = -oo 

g{x,y) = -VG{x,y), (45) 

where G{x,y), the kicking potential, is a determinis- 
tic function of {x, y) which is periodic and sufficiently 
smooth (e.g., analytic), and T is the kicking period. 

The numerical experiments reported hereafter were 
made with the kicking potential G{x,y) = Gi{x)Gi{y), 
where Gi{q) is given by 

Gi (q) = 7: sin 3q + cos q, (46) 
o 

and the kicking period, T ~ 10^®. The number of 
collocation points chosen for our simulations is gener- 
ally TVj, = 10^. In Fig. 4 we plot the PDF of ip ac- 
cording to Eq. (26), and compare it with the numer- 
ical results. In Fig. 5, the moments of ip increments, 
{\p{x2 + Sx) — p{xi)\^), arc calculated numerically as a 
bmction of a; = — a;i| for several values of b (with 
< 6 < 1, and 1 < 6) and its scaling exponents ^f, for 
a; << 1 are checked. The results indicate that with good 
precision (|(5(y9|'') scales with x with an exponent 1 for 
& > 1; otherwise, it scales with x with exponents = b. 
Values of sltc given in Fig. 6. 

The bi-fractal behavior of the exponents is a conse- 
quence of the presence of the domain walls. Indeed, the 
structure function, 

Gb = {\pix2) - pix^)\') , (47) 

for 6 > behaves, for small Ax = \x2 — xi\ as, 

Gb ^ At\Ax\'' + A',\Ax\ (48) 

where the first term is due to the regular (smooth) 
parts of the order parameter p, while the second one 



is contributed by the 0(|Aa;|) probability to have a do- 
main waU somewhere in an interval of length |Aa;|. For 
< 5 < 1 the first term dominates as jAxj 0, while, 
for 6 > 1 it is the second term that does so. 



VI. SUMMARY 

We studied the domain wall-type solutions in the f'^- 
theory in the strong-coupling limit, k' 0, in which 
the equation develops singularities. The scaling behavior 
of the moments of differences of ip, Sip = ip{x2) — ip{xi), 
and the PDF of (p, i.e., Pif), were all determined. It was 
shown that in the stationary state, where the singularities 
are fully developed, the relaxation term in the strong- 
coupling limit leads to an unclosed term in the equation 
for the PDF. However, we showed that the unclosed term 
can be omitted in the strong-coupling limit. We proved 
that to leading order, when \x2 — xi\ is small, fluctuation 
of the ip field is intermittent for 6 > 1. The intermittency 
implies that, Cb = {\(p{xi) — ip{x2)\'' > scales as \xi — 
X2|^'', where ^f, is a constant. It was shown, numerically, 
that for the space scale |a;2 — a^il and 6 > 1, the exponents 

are equal to 1. 



VII. APPENDIX I 

In this appendix we provide a proof of Novikov's theo- 
rem. Consider the general stochastic differential equation 
with the following form, 



d 1 

-g^f = -2L[ip{x, t)] + 77(a;, t) 



(49) 



where L is an operator acting on ip, and 77 is a Gaussian 
noise with the correlation, 

{ri{x, t)Tj{x', t')) = k{x - x')5{t ~ t') . (50) 

The PDF of the random noise has the following form, 

iMv)] = [dv] X 



exp[--J <Fx(fx' dtdt' r]{x,t)B{x - x')5{t - t')ri{x' ,t')] , 
where B{x — x') is the inverse of fc(a; — x'), so that, 

k{x - x')B{x' - x")d'^x = S{x - x") . (51) 



we write the average of ri{x, t)F(r]) over the noise real- 
ization as: 



{7^{x,t)F{i^)) = / 77(x,t)F(77)[dp(77)] 



(52) 



By integrating by parts and using Eq. (51), one finds 
that, 



{7^{x,t)F{i^)) 



d''x"dt"{r-i{x,t)r]ix",t")){ 



dF 



dTjix",t")' 



(53) 



Now, let us assume the function F to have the following 
form, 



^[77] = exp{—iXip{x',t)) 
so that one finds. 



dF 



dipix\t) 



dTj{x",t") d7j{x",t") 
Integrating Eq. (49) with respect to t, we find that 
ip{x\t) = 



(54) 



(55) 



1 



^{x'M)-- dt"L[ipix',t")]+ dt"rj{x',t") . (56) 



ta 

This allows us to show that, 

dv{x\t) ^ 
dr]{x",t") 
1 ^^,„dL[^{x,t"')] 



to 



dri{x' 



L-jl + S{x'^x")0it-t") 



(57) 



where, in the limit t" t, the first term of the right-hand 
side of Eq. (57) will vanish, and we can write 

{T]{x,t) exp[-iXip{x',t)]) 

= {-~i\)k{x - a-')(cxp(-jA(^(a;', t))) . (58) 
where we used, 6{Q) ~ 1 

VIII. APPENDIX II 

In this Appendix we prove that, for example in Eq. 
(24), the relaxation term k'V'^ip makes no contribution 
or anomaly to the PDF of the increments, in the limit 
k' 0. 

The joint probability distribution P{ipi,ip2) satisfies 
the following equation. 



Pt{ipu'P2) = -[(</?! - fl)P]vi 
„3 



[{ip2'ipi)PU+k{0){P^, 



-[-2/c(x)_P(^^ 

, f dXi d\2 , , 

~ / -^-^-^1 exp(iAi(^i -f zA2(;52) 

x{V'^ip{xi) exp[-i\iip{xi,t) - i\2ip{x2,t)]) 

, f dXi d\2 . , . . « 

~ / ^^^2 exp(iAi(^i zA2(^2) 

x(V^(/3(2:2)exp[-iAi(/7(a;i,t) - iA2</3(x2,t)]) . (59) 

The last two terms in Eq. (47) are not closed with respect 
to the PDF. Let us then compute the contribution of the 
unclosed terms. They can be written as, 



©- 

V , 



b = 0.1 
b = 0.5 
b = 1 
b = 2 
b = 3 
b = 4 
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FIG. 5. The moments, {\(p{x2) — (p{xi)\^), as a function of 
\x2 — xi\, obtained via numerical simulation. 
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FIG. 6. Scaling exponents ^6, checked numerically. The 
results indicate that, with a good precision, {|5</5|'') scales with 
X with an exponent 1 for 6 > 1 larger than one, and elsewhere 
scales with x with an exponent ^t, = b. 



= -k' lim V / dxdv]i_^[ipx,x, + Vrf.iV?) 

x5[lp, - ip{x^,t)\ (64) 

where Q.j is the space close to the domain walls. There- 
fore, Eq. (52), in the limit, k' 0, is written as, 

= -k' lim ^'V] [ dxidvli_j^ipxiXiS[fi~'f{xi,t)]. 

3 ^ 

(65) 

Changing the variables from Xi to Zi and integrating over 
dv^d-i^ one finds, 

= -k' hrn^ ^-^^ f edzi^(p^.;,^5[ipi - Lp{zi,t)] 
= lim -^-i-V" / dziip;,.^^S[ipi- (p{z^,t)] 



where e = {2k'Y/'^. Assuming statistical homogeneity, 
we have 



k' NVd- 



= lim 



e y— oo V 



dzi(p^^^J[ipi - (p{zi,t)] (66) 



where N is number of singular lines. Moreover, ^ = 
(fc')^/^, and, ^^'^ is the density of the singular lines 
which, in the limit, V — > oo, is simply the singularity 
density p. Therefore, 



+ OC 



-{k'Y^'^p dz^ifiz^ZiS[(pi - (p{zi,t)] . (67) 



~ik J ^^^^\iexp{i\iipi+iX2'P2) 
X < vV(a;»)e~*^^'^^'^i'*^"'^='^('=^'*) > 

= "^'^ ( / TT^TT^ ea;p(iAi</Ji + iX2'P2) 
J Ztt Ztt 

X < vV(a;^)e"'^''^^''''*^"*^"^^''"*^ >)^, 

= -k'{V^(p{xi)d{(pi - (p{xi,t)) 

xS{lp2 - ip{x2,t)))^. 



(60) 



(61) 



(62) 



Consider one of the terms in the above equation, for 
example, -k' {V'^ip{xi)d[ipi-ip{xi,t)]xS[ip2-'p{x2,t)])^^. 
Assuming ergodicity, it is written as, 

fc' lim / dx,dvli_^\/'^(p{x,)6[ipi -- (p{xi,t)] (63) 

V^oo V Jy 

in the limit, k' limit, only at the points where we 
have singularity this term is not zero. Therefore, we re- 
strict ourselves to the space near the domain walls, 



In the same way in, for example Eq. (21), by changing 
the integration variable from z,; to ipi^ we calculate the 
integral exactly, 

= ^Fj2^{A + 2V2(f>^S{(f), - 1)} 

where A = 27261(1 - 0f )6'(0f - 1) - V26l(l - </)f ). There- 
fore, in the limit, fc' ^ the master equation will give 
Eq. (35). 
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